Ghost propagator and the Coulomb form factor from the lattice 

G. Burgio, M. Quandt, and H. Reinhardt 
Institut fiir Theoretische Physik, Auf der Morgenstelle I4, 72076 Tubingen, Germany 

Abstract 

We calculate the Coulomb ghost propagator G(|p|) and the static Coulomb potential Vbdrl) for 
SU (2) Yang- Mills theory on the lattice. In view of possible scaling violations related to deviations 
from the Hamiltonian limit we use anisotropic lattices to improve the temporal resolution. We find 
that the ghost propagator is infrared enhanced with an exponent Kgh > 0.5 while the Coulomb 
potential exhibits a string tension larger than the Wilson string tension, ac — 2(7. This agrees 
with the Coulomb "scaling" scenario derived from the Gribov-Zwanziger confinement mechanism. 
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I. INTRODUCTION 



Yang-Mills theories in Coulomb gauge have attracted increasing interest in the last years, 
both in the continuum and on the lattice. Non-perturbative analytic predictions based on 
the Gribov-Zwanziger confinement mechanism [H [2] turn out to be concise and elegant in 
this gauge and numerical simulations have confirmed them so far. 

In covariant gauges the Gribov-Zwanziger approach requires a restriction of the functional 
integral beyond the Faddeev-Popov mechanism [S], which has long been known to be insuffi- 
cient to extract a unique field representative along the gauge orbit and thus cannot be used 
to define the partition function beyond perturbation theory [H HJ [5]. As a remedy one adds 
extra terms to the action which are, in general, non-local and often referred to as the Hori- 
zon function or Horizon condition. The purpose of these extra terms is to limit the fields 
in the functional integral to the so-called fundamental modular region A, which contains 
only absolute minima of the gauge fixing functional and thus eliminates the over-counting 
of gauge copies from the same orbit (except for topologically non-trivial copies). 

In practice, however, it is almost impossible to limit the functional integral beyond the 
so-called first Gribov region f2 D A, where the Faddeev-Popov operator is positive defi- 
nite. Moreover, any restriction on the integration range imposed other than through the 
Faddeev-Popov mechanism will break the Becchi-Rouet-Stora-Tyutin (BRST) symmetry 
[6l [7] . Whether such enlarged action including Horizon terms exhibits some other symmetry 
and what the possible consequences in covariant gauges are (e.g. in terms of contributions 
from non-standard condensates [HHH]) is still under active debate, cfr. Ref. [121 [13] and 
references therein for recent results. 

In contrast to the situation in covariant gauges sketched above, the physical implications 
of the Gribov-Zwanziger approach are quite transparent in Coulomb gauge. In the Hamil- 
tonian formulation p^HTT] . for instance, the Gribov-Zwanziger idea amounts to a mere 
projection of the physical Hilbert space onto the subspace of states satisfying GauB's law 
[T8| [19] . Besides being conceptually simpler, this has at least two advantages: 

• Unlike the Kugo-Ojima approach [20] in covariant gauges, where the existence of a 
globally conserved BRST charge QBRST|^)ph = is essential, the Hamiltonian ap- 
proach in Coulomb gauge requires no additional assumptions to ensure GauB's law, 
i.e. a vanishing colour charge on physical states, (5c|^)ph = 0. 
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• The construction of the physical Hilbert space in Coulomb gauge is much simpler than 
in the fully gauge-invariant Hamiltonian approach [21]. In fact, the consequences of 
GauB's law in Coulomb gauge can be implemented exactly in a functional integral, 
which is therefore well suited for further approximations without the impediment of 
additional constraints such as the conservation of Qbrst- 

The Horizon condition in Coulomb gauge implies both a static (i.e. equal-time) transverse 
gluon propagator [1] which vanishes at zero momentum, and a ghost dressing function which 
diverges in the same limit [2]. Physically, the last statement means that the Yang- Mills 
vacuum behaves as a perfect color dia-electric medium (i.e. a dual superconductor), because 
the dielectric function of the Yang-Mills vacuum agrees with the inverse ghost dressing 
function p2J. As a consequence the dual Meissner effect, which has long been a model 
for the origin of the confining force in so-called Abelian gauges, can also be applied to 
Coulomb gauge. A natural quantity to study confinement in Coulomb gauge is the non- 
Abelian Coulomb potential, which provides an upper bound for the quark- ant iquark free 
energy V{r) < Vc{r), i.e. there is no confinement without Coulomb confinement |23j . 

For all these reasons Coulomb gauge is best suited for direct investigations of the QCD 
wave functional. After pioneering analytical and numerical work [27], such studies 

have experienced a broad popularity in the literature, cfr. [28H37j. In particular, the Hamil- 
tonian approach lends itself to variational formulations [I5l [T71 EI] , which allow to address 
non-perturbative problems in Yang-Mills theory in a rather direct and concise way. A main 
ingredient in these techniques are static (equal-time) two-point functions, so that a direct 
investigation of the Gribov-Zwanziger scenario for such correlators is important. 

There have been a number of studies of this subject in the literature |38lflT] . From 
previous analyses of the gluon propagator in Coulomb gauge [i2] - H5] it is known that the 
difficulties in reaching the Hamiltonian limit on a lattice with a finite time-resolution can 
have drastic consequences: the resulting scaling violations, if untreated, prevent a multi- 
plicative renormalization of the gluon propagator and also affect the correct analysis in the 
deep infrared. We therefore compute the Green functions on anisotropic lattices with a high 
temporal resolution and study the effects of approaching the Hamiltonian limit. We also 
employ a high quality of gauge fixing, particularly important for the ghost propagator. 

The paper is organized as follows: in the next section, we discuss the exact definition of 
the relevant correlators, both in the continuum and on the lattice. We will also display the 
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expectations on the conformal behaviour of the correlators in the deep-infrared, as raised 
by variational and functional methods. Section 3 contains a description of our numerical 
setup, as well as a detailed discussion of our results for both the two-point functions and the 
Coulomb potential. We also compare to continuum methods and give arguments to explain 
possible deviations. Finally, we conclude in section 4 with a brief summary and outlook. 



II. CORRELATORS IN COULOMB GAUGE 
A. Correlators in the continuum 

At fixed time t, the Coulomb gauge condition V ■ A = is complete (up to possible 
singularities and Gribov ambiguities, see Ref. i-e. any residual gauge symmetry is 

space-independent and thus acts as a global gauge at fixed time t. In this situation, we are 
interested in the static (i.e. equal time) propagators, 

^'p> = Yw: ''-"''-p- *)> = iilf / 1^ {Anp)A%-p)) (1) 

G(P) = = r i^iPy'i-P)) = tr ( (-D ■ V)-' > (2) 

IPI ^^A 

^o(p) = ^(^S(P,^K(-P>^)> ~ (3) 
~ Vciv) = / tr ((-D ■ V)-^(-V^)(-D ■ V)-^) , (4) 

where in the last equation contributions from non contact terms have been dropped [47J . A^^ 
is the dimension of the adjoint colour representation and tr is the normalized colour trace; 
we use the colour group G = SU (2) with A^^ = 3 throughout this paper. Moreover, the 
covariant derivative in the adjoint representation reads D"^ = V 5^*^ -|- ge°'^^A'^, and (— D ■ V) 
is the Faddeev-Popov operator in Coulomb gauge. 

Within the Hamiltonian variational approach one finds two sets of solutions for the gap 
equations, called "critical" [El HE] and "subcritical" |19]. Here we will be mainly interested 
in the critical solution, which is also believed to be the physically relevant one. In this case 
one finds a Coulomb potential Eq. (|4]) which rises linearly as a function of the distance |18] . 
Furthermore, the gluon and ghost propagators in Eq. 
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Tp) exhibit a conformal scaling in 



the deep infrared, i.e. they behave as a power of momentum |p| ^ 1, 

^(P) ~ Ar: , dip) ~ ■ (5) 

ppgl P gh 

Similarly, asymptotic freedom indicates that the large momentum behaviour of the prop- 
agators essentially follows from the bare action, modified by logarithmic corrections with 
appropriate anomalous dimensions, 

""^^^ ^ |p|log-..M ' ^^P) ^ ■ 

In both regimes, the relevant exponents are further constrained by the assumption that the 
(static) ghost-gluon vertex is essentially trivial, i.e. proportional to the bare vertex, for any 
kinematical configuration. This leads directly to the so-called sum rules 



7gi + 27gh = 1 . (7) 

The two solutions of the variational approach are derived under the same assumption, with 
the specific values for the exponents [T71 H81 - I50] : 



Kgi - 1 


/tgh — 1 


(8) 


K^i ~ -0.6 


~ 0.8 


(9) 


7gi = 


7gh = 2- 


(10) 



Notice that Eq. ^ implies, in general, an infrared mass generation for both gluon and ghost, 
unless Kgi = 1 and /tgh = 0, their tree- level values. Such an infrared behaviour agrees with the 
original analysis of Gribov [Ij and, in particular, implies an infrared vanishing static gluon 
propagator D{p) and an infrared divergent ghost form factor d{p) — )■ oo, as |p| 0. 
Physically, these findings translate into a diverging gluon self-energy uj{p) ~ D(p)^\ and a 
vanishing dielectric function of the Yang-Mills vacuum, e(p) ~ d{p)^^, which in turn implies 
dual superconductivity |22] . In Ref . the ghost and gluon propagator were also calculated 
from the renormalization group flow equations, assuming a bare ghost-gluon vertex and UV 
propagators with vanishing anomalous dimensions (7gi = 7gh = 0); the corresponding IR 



exponents Kgi ~ —0.28 and Kgh — 0.64 turn out to be somewhat smaller than in Eq. (10). 

The large momentum behaviour cited above for both propagators shows interesting non- 
perturbative effects as well. Despite asymptotic freedom (and in contrast to Landau gauge), 



the predictions Eq. (10) for the anomalous dimensions cannot be directly compared to per- 
turbation theory [321 ES], since no obvious rainbow-ladder-like re-summation technique is 
available in Coulomb gauge. In fact, it is not even clear whether static propagators can be 
accessed in standard perturbation theory in the first place: the Slavnov- Taylor identities im- 
ply a highly non-trivial energy (i.e. po) dependence in the Green's functions [5lj, which could 
spoil the naive static limit and thus require a non-perturbative treatment anyhow [551 - 157] . 
Alternatively, a perturbative approach based on the functional renormalization group has 
been attempted in Ref. [HHlEn]; their prediction 7gi = 3/11 and 7gh = 4/11 differs, however. 



quite substantially from the non-perturbative results in Eq. (10). 



B. Correlators on the Lattice 

A comparison between continuum and lattice Coulomb propagators was pioneered in 
Ref. [6^; later studies with different techniques gave, however, contradicting results [38lflT] . 
The difficulty lies in the fact that most continuum results for static quantities are naturally 
obtained in the gauge Aq = 0, which cannot be attained directly in lattice calculations 
because of the compactified time direction. Furthermore, any finite lattice has a finite time 
resolution and the associated discretization artifacts preclude the Hamiltonian limit and give 
rise to severe renormalization problems. A way to circumvent these issues was first proposed 
in Ref. [^2] - H^ . There it was shown that the static gluon propag ator of Eq. Q agrees, after 



dealing with compactification and renormalization artifacts, with the exponents in Eq. (10). 
More precisely, the static gluon propagator in = 3 + 1 was shown to be surprisingly well 
described, over the whole momentum range, by Gribov's original proposal [T]: 

D{p) = , . (11) 

As argued in Ref. [121133], the problem in the "naive" extraction of -D(p) from the lattice data 
lies in the standard, isotropic discretization of Yang-Mills theories. The use of anisotropic 
lattices, which approach the Hamiltonian limit pTj for large anisotropy ^ = ag/at (see 



Sec. Ill A), was therefore proposed as a general tool in Coulomb gauge investigations. Fur- 
ther independent studies in Ref. [621 IBS] have confirmed the presence of such effects and the 
improvement of lattice results for ^ > 1 for the gluon correlators. Incidentally, the authors 
there were still not able to extract a coulomb string tension from the analysis of the temporal 



propagator, Eq. ([s]); the supposed equivalence between Eq. (|3| and Eq. ^ [21 [23] remains 
therefore an open issue. 

In this paper we will continue and complete the lattice analysis of Green's functions in 
Coulomb gauge for pure SU (2) Yang-Mills theory at T = 0, which was started in Ref. jlJUH]. 
We will concentrate on the ghost form factor d{p) and the Coulomb potential Vc'(p), and 
study, in particular, the anisotropy effects and the validity of the confinement scenario 
sketched above. A study of the Coulomb gauge quark propagator with light dynamical 
fermions has been published elsewhere |45j . 



III. RESULTS 



A. Numerical setup 

Pure Yang-Mills theory, as any continuum field theory, can be formulated in the Hamilto- 
nian picture [6^, where space and time are treated separately. Upon discretization, simple 
RG-group arguments [65] indicate that for any given isotropic lattice version of the field 
theory in question a corresponding anisotropic counterpart lying in the same universality 
class exists [66j. The latter describes the same physics, albeit with two different lattice 
spacings as and at for the space and time directions. For Yang-Mills actions in D = d + 1 
dimensions built in terms of the character x of m x n Wilson loops Pfj,u{x; m, n), 

Six; m,n) = P^J2J2 diioO^ ^)]) (12) 

the anisotropic counterpart will read (see Ref. [61] and references therein): 

S'{x;m,n) = J2^(3s (l ~ dim(x) ^^^''^^''"'''^0 

For each choice of (3^ ^ (if the two lattice spacings and at have to be determined non- 
perturbatively. The couplings are usually parametrized as = ■ 7 and [3f = (5^ ■ 7"^, 
where is a common coupling factor, while 7 is the hare anisotropy; it is related to the 
true (renormalized) anisotropy = as/ at through the renormalization constant 77, i.e. we 



have ^ = 'y ■ rj. The non-perturbative value of t] can be shown to slowly approach, in 
the weak coupling limit, its perturbative expression [UTl [U7] . In much the same fashion, the 



transition from an isotropic action such as eq. (12) to its anisotropic counterpart eq. (13) can 
be generalized to lattice actions of the form S = J2xmn^x,mn S{Xj''T^y^)- Such anisotropic 
models are very useful in lattice applications, see e.g. Ref. [6ll |68] and references therein. 

Turning to our specific case, we will concentrate on the standard Wilson one-plaquette 
action in 3+1 dimensions, which for SU{Nc) pure gauge theory reads: 

Sw = J2^ 11 i^'W^^ iP,Ax))]) , (14) 
while its anisotropic version is given by: 



Ux))] (15) 



Sw = J2\PsJ2 JT^' [Tr (P.,(x))]) + A (l - ^Re [Tr (P. 

X I i>j=l V c / \ c 

{3 3 
^l(^-jr^^ (^.. W)]) + ^ E (l - [Tr (^.4(x))] 

i>j=l ^ / ' j=l \ 

The values for at and = atC, for various choices of /3 and 7 have been extensively studied 
in the case G = SU{3), see e.g. Ref. |63|, |69]. For SU{2) we are only aware of one study 
from the literature [7Q|. In both cases non-perturbative effects turn out to be quite large, 
so that the analytic one-loop calculations for r], which are in principle available for any A^^, 
(3 and ^ [6IJ, cannot be trusted for practical applications [671 169]. 

Since we will concentrate on the case of two colours, = 2, we have decided to re- 
determine the relevant parameters independently by imposing rotational invariance for the 
static potential extracted from space- and time-like Wilson loops [HZl EH] • Our best estimates 
for 7 and are given in Tables |T] and [IT] for a selection of values for /3 and ^. These results 
have been obtained on x (^ L) lattices of spacial extension up to L = 32 and anisotropies 
up to = 4; we have sampled ^(1000) configurations using a heath-bath algorithm combined 
with over-relaxation, which for SU{2) can be easily extended to the anisotropic case. The 
same algorithm was also used to generate the configurations on which the measurement of 
Green's functions described in the next section have been performed. 

The results for 7 are mostly compatible within errors with those of Ref. [70], while we 
find some discrepancies in the scale a^. We have checked that our predictions for = at 
in the isotropic case C = 1 agree with others in the literature (cfr. e.g. Ref. [71]), while 
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the results of Ref. [70] are always higher than ours at the lower end of the scaling window 
f3 < 2.3. Such discrepancies could be due to different systematics inherent to the method 
used; higher precision might be needed to settle the issue. 



/3 


2.15 


2.2 


2.3 


2.4 


2.5 


2.6 


2.7 


e = 2 


1.654(3) 


1.672(3) 


1.712(3) 


1.754(4) 


1.796(5) 


1.835(6) 


1.870(9) 


e = 3 


2.375(3) 


2.407(3) 


2.474(4) 


2.545(4) 


2.608(5) 


2.663(7) 


2.710(9) 


e = 4 


3.106(4) 


3.151(4) 


3.243(5) 


3.331(5) 


3.406(6) 


3.466(7) 


3.511(9) 



TABLE I. Bare anisotropy 7 vs. true anisotropy ^, for various couplings f3 in SU{2). 





2.15 


2.2 


2.3 


2.4 


2.5 


2.6 


2.7 


e = i 


1.196(6) 


1.061(6) 


0.821(6) 


0.616(6) 


0.441(5) 


0.290(5) 


0.231(5) 


e = 2 


1.355(7) 


1.206(6) 


0.940(6) 


0.711(6) 


0.511(5) 


0.338(5) 


0.268(5) 


e = 3 


1.391(7) 


1.239(7) 


0.968(6) 


0.732(6) 


0.526(5) 


0.345(5) 


0.270(5) 


e = 4 


1.406(8) 


1.254(7) 


0.979(6) 


0.739(6) 


0.530(5) 


0.346(5) 


0.271(5) 



TABLE II. Spatial lattice spacing for the choices of 7 as in Table |I| Data are in GeV , assuming 
a Wilson string tension of o" = (440 MeV)^. The first line gives our best estimate for the isotropic 
scale, cfr. e.g. Ref. [7T] . 

It should be noticed that the scale raises with ^ at a fixed coupling /3, up to ^ > 3 
where a plateau is reached. In order to simulate at the same physical point for observables 
involving spatial links only (like the static correlators we are interested in), one thus needs 
to tune (3 to C, as the latter increases. 

B. Coulomb gauge and the lattice Hamiltonian limit 

In Ref. [12], HS] , it was shown how to quantify and treat the explicit at dependence appear- 
ing in static observables after integrating their non-static counterparts w.r.t the energy pq. 
However, there are still more subtle effects of the finite time resolution in lattice simulations. 
In Ref. [68] it was shown that the lattice Yang-Mills spectrum suffers from discretization 
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effects of order ^{a^) up to i^(aj), depending on the quantum numbers at hand. Such effects 
are purely dynamical and linked to the rate at which the theory on an anisotropic lattice 
reaches the physical Hamiltonian limit ^ — )■ oo in which the eigenstates and the spectrum 
are well defined. Lattices with ^ = 5 were eventually used in Ref. [BBJ to extract the masses 
of the physical states. 

The goal of our investigation is to compare Coulomb gauge lattice results with continuum 
Hamiltonian calculations. Discretization effects at least of the same order as for the glueballs, 
if not larger, should therefore not come as a surprise. Let us illustrate this effect using the 
Coulomb gauge functional as an example. 

To fix Coulomb gauge on each MC generated configuration we adapt the algorithm of 
Ref. [721 [73] as described in Ref. [i2tili] . i.e. we first fix one of the = 2^ = 8 center 
flip sectors and then maximize, via a simulated annealing plus ensuing over-relaxation, 
separately at each time slice t the gauge functional 

Fit) = f^f (x, t) , t/f (x, t) = gi^, t) f/,(x, t) + 1 1) (16) 

with respect to local gauge transformations g{'^,t) G SU{2). To further improve the gauge 
fixing quality, we start the gauge fixing engine with random gauge copies of the initial 
time slice (with usually 5 < -n-c < 10), and select the gauge- fixed copy with the highest 
value of F(t). We then combine these best gauge-fixed time slices into a best gauge- fixed 
configuration for the flip sector chosen. Finally, we proceed to the next sector and repeat 
the procedure, selecting in the end the sector which gave the highest global functional 
F = '^fF{t). Every time slice in the final gauge-fixed copy for each configuration has 
therefore been selected among at least Uc ■ = 40 ... 80 gauge fixing runs with random 
starting points from all topological sectors. On top of this elaborated procedure to fix 
the spatial gauge freedom in each time slice separately we fix the residual temporal gauge 
symmetry g(t) via the integrated Polyakov gauge, as described in Ref. [i2] - H^ . 

As we have seen in the previous section, a constant spatial cut-off for increasing 
anisotropy ^ can be obtained by tuning /3 to ^; = a,, will then decrease linearly with 
Naively one would then expect all observable which only depend on the spatial links 



to be independent of C,. For instance, the Coulomb gauge functional Eq. (16) should, in 
principle, be a function of only, at least if our algorithm is good enough to approach 
the absolute minimum reliably, and if finite volume effects can be ignored [71]. This is not 
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the case, as can be seen from Figure [T] the left chart displays the best value of the gauge 
functional F from our algorithm as a function of the anisotropy for three fixed values of 
the lattice spacing = 1.060(6), 0.556(5) and 0.350(5) GeV~^. As the Hamiltonian limit 



a 0.87 



0.83 
0.5 



H 1 

a^ = 1.060(6) GeV 




g Q a^ = 0.556(5) GeV-' 




□ = 0350(5) GeV-1 - 




a 




H 

B 

e 


Li- 


e 

e 


Li."^ 




10"^ 



2.5 



/ ^ 


/ G 




' . \^ 1.060(6) GeV"'' 




^ a^- 0.556(5] GeV"^ 




^ = 0.350(5) GeV"^ 



(a) (b) 

FIG. 1. (a): dependence of the gauge fixmg functional F^{a^ on the anisotropy ^ at fixed spatial 
cut-off ag. (b): Deviation of the gauge fixing functional from the Hamiltonian limit, F^{a^ — 
Fooio-s), as a function of the temporal lattice spacing at, together with its leading power corrections. 



is approached by increasing ^, the functional F decreases, i.e. the configurations at fixed 
time slice become "rougher" although ag is kept constant; as can be seen from the figure, 
the corrections with ^ are several orders of magnitude higher than the Gribov noise. The 
leading order in the corrections to F^ are well described by a power law in oc at- Our 
best estimates for the asymptotic values F^o together with the coefficient of the leading order 
corrections are given in Tab. Ill, while the right chart of Fig. [T] shows F^ — F^o versus the 
leading correction in at- For the stronger couplings, this leading correction is ff{af), while 
for weaker couplings we find it to scale as ^{af). 

This sensitivity to ^ (or at) does not only occur in the gauge fixing functional. In 
Ref. IISI ESI E2j a strong dependence on the anisotropy was also seen for the temporal 
gluon propagator Eq. ([s]), although this was numerically shown to be energy (i.e. Uq) inde- 



pendent |l2j. Also, preliminary results indicate that the Gribov mass M in Eq. (11) slightly 

we 



increases with ^, saturating to a value M ~ 1.2(2) for ^ = 3 — 4. In Sec. IIIC and HID 



will see that the ghost propagator Eq. ^ and the Coulomb potential Eq. (|4]), which are by 
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7? 

^ oo 


C2 


C4 


as = 1.060(6) GeV-^ 


0.8333(2) 




0.152(1) 


as = 0.556(5) GeV-^ 


0.8620(2) 


0.058(1) 




as = 0.350(5) GeV-i 


0.8737(2) 


0.055(1) 





TABLE III. Asymptotic extrapolation of together with the coefficient of the leading power 
correction of oder ^(a^) (a^ = 1.060(6) GeV-^) and ^(a^) (a, = 0.556(5) and 0.350(5) GeV-^). 

definition independent of the time-like links, both turn out to be sensitive to the anisotropy. 

From Fig. [l] one can infer that lattices with ^ > 6 would be needed to minimize the 
at corrections in F, i.e. to reach the plateau at F^o- Since we also need a large spatial 
extension L ~ 50 to explore the deep infrared momentum region, the combination of these 
two requirements would force us to simulate on lattices of temporal extension Nt ~ 300, 
beyond the computational power at our disposal. Moreover, not all observables need to 
be as sensitive to ^ as the gauge functional. We will therefore restrict our investigation to 
^ = 1 ... 4 and attempt to extrapolate to larger ^ wherever necessary. Of course, a direct 
confirmation of our results on large lattices with high anisotropics would be highly desirable. 

C. Ghost Form Factor 

The ghost form factor in Coulomb gauge, Eq. ([2]), has been discussed in Refs. [SS j I1T | [7Sl 
[75] ; neither its ultraviolet nor its infrared behavior could be determined conclusively. In the 
UV, the primary goal is to check the sum rule for the anomalous dimensions in Eq. ([T]). In 
the IR, the main question is whether the ghost propagator is compatible with an infrared 
finite behaviour, as is the case in Landau gauge |77j. If this were true, it would of course 
spoil the Gribov-Zwanziger mechanism and the dual superconductor argument of Ref . [22] . 

We calculate d{p) by inverting the Coulomb gauge Faddeev-Popov (FP) operator through 
a conjugate gradient algorithm on lattices of spatial sizes up to 54^ and anisotropics up to 
^ = 4. Although for most configurations the algorithm works quite well, there are some 
"exceptional" cases where very small eigenvalues of the FP operator make the inversion 
ill-conditioned, signaled by very bad convergence. On the other hand, we expect exactly 
these configurations to contribute substantially to the infrared divergence of the ghost form 
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factor, since they lie close to the Gribov horizon. Our current procedure is to exclude such 
configurations from the Monte-Carlo ensemble; this obviously creates a bias in the data 
which potentially suppresses the ghost propagator at very low momenta. 

Whether or not these near-horizon configurations have a measurable effect on the ghost 
propagator also depends on the frequency with which they appear within a MC sequence; 
contrary to Landau gauge [7^, this seems to be relatively stable upon improvement of the 
gauge fixing. In our studies we have observed an ill-defined FP operator in about one of 
every 300 configurations, with very large uncertainties (i.e. there were also MC runs with 
over 1000 configurations that did not show a single near-horizon configuration). Preliminary 
studies of individual singular configurations with improved pre-conditioners and numerics 
indicate that their contribution to the ghost form factor can be enhanced by up to a factor 
40 as compared to the ensemble average, again with large uncertainties due to the bad 
condition number. Clearly, this is a zero measure times infinite contribution problem that 
can only be decided by much improved statistics combined with better inversion algorithms. 
In this study we have not been able to resolve this issue quantitatively; the exponent which 
we will extract in the following must therefore be considered as a lower bound for the correct 
infrared behaviour. 

Results for d{p) at ^ = 1 are shown in Fig. |2](a), while in Fig. we give the coeffi- 
cients Z{I3) needed to multiplicatively renormalize the data at different coupling, arbitrarily 
scaled such to let d{p) ^ 1 at large momentum. Although at first glance the data agree 
qualitatively with the above expectations, on a closer look two problems appear. First, 
small scaling violations can be measured within errors in the ultra-violet region; a fit to a 
logarithmic behavior is thus afflicted with large errors. Second, the infrared region shows 
two different behaviors for an intermediate momentum range 0.5 GeV< |p| < 1.2 GeV and a 
low momentum range 0.2 GeV < |p| < 0.5 GeV. In the first regime the data are compatible 
with a power-law behavior 

d{\p\) oc |p|-^ (17) 

with K, ~ 0.5, while in the low-momentum region the effective k decreases, approaching 
a value k, < 0.4. By going to higher anisotropics this behavior softens and the data be- 
come more and more compatible with a stable power-law. To illustrate this effect, in 
Fig. |3](a) we highlight the infrared behavior showing |p|'^™c?(p) for different anisotropics, 
where Km = 0.373(6) is the effective exponent we have measured for L = 54, ,^ = 1 in the 
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(a) Cb) 

FIG. 2. (a): Ghost form factor (i(|p|) for anisotropy = 1. (b): Renormalization coefficients for 
tfie gliost propagator, as a function of the coupling (3. 

lowest momentum region. The data for ^ = 1 go to a constant while for higher anisotropies 
a power-law still describes the data well. In Fig. [sj^b) we show the IR exponents n obtained 
from three different fitting methods for each set of data at fixed anisotropy ^, i.e. tempo- 
ral cut-off af. the lower curve corresponds to a pure Ip]^" behaviour cutting the data at 
IpI < 1.5 GeV; the middle one adds a linear subleading correction and leaves the same 
cut; the higher one combines the linear correction to the power law with a cut at |p| < 1 
GeV. Fitting all these values of k to a constant Kgh, constrained to be the same for all three 
methods used, plus power law corrections in at we find 



Kgh = 0.55(1) 



(18) 



with leading correction of order ^(a^) and a x^/d.o.f. ^ 0.3. 

Because of the cut on exceptional configurations discussed above we have of course in- 



troduced a bias in our data. Therefore the value in Eq. (18) can only be considered as a 



lower bound on the actual value of Kgh. Still, as far as we know, this is the first time that 
an IR enhanced ghost form factor can be reliably proven to exist on the lattice, making the 
Gribov-Zwanziger confinement scenario in Coulomb gauge fully consistent. Indeed, while 
in Ref. |S] it was shown that the Gribov gluon propagator Eq. (11) is IR equivalent to 



the massive Landau gauge solution, in the latter case lattice simulations always give an 
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■ 0.5 1 1.5 2 2.5 3 3.5 . « .. ^ (GeV)"'' 

iPl (GeV) 

(a) (b) 

FIG. 3. (a) Infrared behavior of |p|'^'" '^(p) ^or different anisotropies ^. (b) k extrapolation for 
at — )• 0; Different symbols correspond to different fitting methods to estimates k; points at fixed at 
(^) correspond to the same data set for ci(|p|); lines give the leading correction oc at to Kg^. 

IR finite, i.e. tree-level like ghost form factor. How this can be made to agree with the 
Gribov-Zwanziger mechanism is still a debated issue. 

In general, although any divergence {ugh > 0) of the ghost form factor would be sufficient 
to support the Gribov-Zwanziger mechanism, the exact value of the infrared exponent Kgh 
matters for the sum rule Eq. ([T]) , which is violated by our findings combined with the results 
for Gluon propagator in Refs. [12H13]. There are several possible resolutions to this issue: 
first, our data for the ghost propagator only covers the range down to |p| ~ 200 MeV 
with sufficient statistics and precision. Further changes in Kgh at a much smaller scale 
\iR <^ Aqcd can not be ruled out, although we currently do not see the onset of such 
deviations, while the appearance of an additional small scale in pure YM theory would raise 
a different kind of interpretation problem. Moreover, as we have explained above, our data 
for is only a lower bound since we have neglected singular configurations from near 
the Gribov horizon. Whether or not these configurations contribute substantially to the 
ghost propagator and its exponent Kgh cannot be predicted with our current computational 
resources. Most likely, it would also require algorithmic changes in the inversion method for 
ill-conditioned FP operators, and maybe even a change in the fundamental MC generator 
to better sample the near-horizon region in field space. 
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Besides these caveats at our numerical data, the origin of the sum rule Eq. ([7]) itself leaves 
room for discussion. In essence, it is based on the assumption that (i) the gluon and ghost 
propagators have conformal (power-like) behaviour in the infrared and (ii) the ghost-gluon 
vertex receives no radiative corrections (other than an overall multiplicative renormalization) 
at low momenta, i.e. it is essentially bare. Both results are borrowed from Landau gauge, 
where Taylor's theorem gives a firm explanation why the ghost-gluon vertex is trivial as 
p — 7- 0. The ensuing sum rule for gluon and ghost propagator are then simple consequences, 
at least as long as massive solutions in the infrared can be ruled out, for which the ghost 
and gluon behaviour would decouple. All these assumptions are well confirmed by lattice 
simulations in Landau gauge. 



For Coulomb gauge, however, the situation may be more involved. A careful analy- 
sis of the Slavnov- Taylor identities in Coulomb gauge exhibit a complicated interplay be- 
tween transversal (spatial) and longitudinal (temporal) degrees of freedom, see Eq. (4.12) 
in Ref. [51]. Integrating over all energies po to reach the equal-time limit could therefore 
induce non-trivial structure in the static Green functions, even though the proof of Taylor's 
theorem carries over to Coulomb gauge in the limit |p| — )■ at any fixed Pq. Unfortunately, 
no analytical calculations of the static vertex at low momenta based on the Slavnov- Taylor 
identities has so far been possible, and the corresponding calculation on the lattice have 
not been carried out. In Ref. [78], within the Hamiltonian approach, the Dyson-Schwinger 
equation for the ghost-gluon vertex was solved at the one loop level and little dressing was 
found. Furthermore, Landau gauge calculations in three (euclidian) dimensions also exhibit 
little dressing of the ghost-gluon vertex; it is however unclear if these results simply carry 
over to the four-dimensional static quantities in Coulomb gauge. 



Let us now consider the ultra-violet behavior, where exceptional configurations play no 
role and the situation is much clearer. As mentioned above, slight scaling violations can 
be measured in the data for ^ = 1. These diminish as ^ is increased and the anomalous 
dimension of the ghost field can be determined in the Hamiltonian limit. In Fig. |4] we 
compare d{p) for the lowest and highest anisotropics ^. For the latter we show our best 



fit to the expected asymptotic behaviour with 7 given by the sum-rule in Eq. (10); the 
ultra-violet mass scale from the fit is m = 0.21(1) GeV, cfr. Eq. ([6]). 
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FIG. 4. UV behavior of d~^{p) for different anisotropics ^ compared with Eq. ([g]). 
D. Coulomb Potential 

The Coulomb potential, Eq. Q, has been intensively investigated in the literature [62|[79l~ 
IM] . While there is general agreement that its infrared behavior is determined by a Coulomb 
string tension larger than the Wilson string tension |23j , the value quoted for the ratio ctc/ct 
varies among the above works. 

We have calculated the Coulomb potential Eq. Q for lattices of spatial extension up 
to 40^ and anisotropies up to ^ = 4. Exceptional configurations did not play a major role 
in this case; this may be either due to the smaller volumes and lower statistics employed, 
since the computation of Eq. Q is much more expensive than Eq. (|2|, or to the better pre- 
conditioning due to the Laplace inversion, or both. In Fig. [5] we show the infrared behavior 
of IpI^ Vc'(p)/(8'7rcr) for different values of ^; the reason for the normalization will be clear in 
the following. As in the case of the ghost form factor, we observe a change in the data in the 
Hamiltonian limit. The corrections due to the finite temporal resolution seem to saturate at 
our highest values of C,, although error bars are still quite large. Also, the bending down of 
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FIG. 5. Infrared behavior of |p|^ Vc'(p)/(87r(T) for different anisotropics ^. 

the data at ^ 0.5 - 1 GeV, observed in all above quoted analyses for ^ = 1, becomes milder, 
although it does not disappear. Its presence should indeed not come as a surprise: at least 
for the physical potential V{r), which we know to be given by |85j 



V{r) 



IT 



1 



we expect, after introducing a suitable IR cutoff A and taking A — 

^2 



\prV{p) = Sna 



TT 



|pr + 



(19) 



(20) 



since the contribution from /i in Eq. (19), which gives a term oc |p| in Eq. (20), will vanish 
with A. Thus |p|^ ^(p) must have a minimum as |p| — ?■ 0; since its asymptotic UV behaviour 
is also, up to logarithmic correction, oc |pp, we only have two possibilities: either |p|^V^(p) 
is a monotonic function of |p| or, if it has some local minimum for some |p|m > 0, then it 
must have a maximum in < |p|Af < |p|m- 

If Vc asymptotically behaves like ac r the data in Fig. |5] must approach a constant as 
|p| — )■ 0, giving direct access to the Coulomb string tension. It is natural to parametrize 
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them as in Eq. (20): 



+ a|p|+7|pr + ^(|p| )• (21) 



Svro" a 

where we would expect 7 > and a — )• in the thermodynamic hmit. Fitting the data 



to Eq. (21) and introducing different cuts at different momenta we obtain for the Coulomb 
string tension 

crc = 2.2(2)a. (22) 

Notice that, even though we don't have much data in the IR region, if we optimistically 
constrain a = we obtain a higher value for the Coulomb string tension, ac = 2.5(1) o". 
Better statistics and data from larger lattices and/or anisotropics in the low momentum 
region would of course be welcome to improve the result. 

IV. CONCLUSIONS 

We have shown that a sufficient approach to the Hamiltonian limit, which on the lattice 
translates into anisotropic lattices with a high temporal resolution, is crucial for a successful 
investigation of the ghost form factor and the Coulomb potential. The effect of the anisotropy 
on static correlators lies in the different dynamics for fixed spatial lattice spacing as at 
decreases, as we have shown in Sec. IIIB[ This explains the dependence on ^ found also for 



explicitly energy independent observables that can be evaluated in a single time slice. 
In Sec. [Ill C we have shown that the infrared exponent of the ghost form factor saturates 



at ~ 0.55(1) when taking the Hamiltonian limit. Although our estimate may only be a 
lower bound due to the contribution of near-horizon configurations, the infrared divergence of 
the ghost form factor confirms the Gribov-Zwanziger [H |2j confinement mechanism and the 
vanishing of the dielectric function for the Yang- Mills vacuum [22] . We have also confirmed 
quantitatively the sum rule for the anomalous dimensions of the static gluon and ghost 
fields in Coulomb gauge. The similar sum rule prediction in the infrared would require 
= 1 in order to be compatible with the Gribov formula confirmed in ref. jSj. This 
is clearly violated by our best fit Kgh ~ 0.55(1). We have discussed possible arguments to 
resolve this discrepancy, by critically analyzing both our data and the origin of the sum rule 
itself. A conclusive settlement of this issue must be left to a future investigation. Finally, in 



Sec. |III D] we have shown that the Coulomb potential is linearly rising in position space, and 
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the corresponding Coulomb string tension settles nicely in the Hamiltonian limit ^ — > oo, 
where we extract a value of ac — 2.2 a. 

As indicated above, future investigations in this direction should concentrate on the ghost 
form factor and Coulomb potential on much larger spatial lattices to probe lower momenta, 
with improved statistics and better numerics to tackle the rare near-horizon configurations. 
This may help to settle the remaining question of the infrared sum rule for the exponents 
in the power- law of the low order Green functions. 
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